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ABSTRACT 

The astrophysical origin of r-process nuclei remains a long-standing mystery. Although some 
astrophysical scenarios show some promise, many uncertainties involved in both the astrophys- 
ical conditions and in the nuclear properties far from the /3-stability have inhibited us from 
understanding the nature of the r-process. The purpose of the present paper is to examine the 
effects of the newly-derived microscopic Hartree-Fock-Bogoliubov (HFB) mass formulas on the 
r-process nucleosynthesis and analyse to what extent a solar-like r-abundancc distribution can 
be obtained. The r-process calculations with the HFB-2 mass formula are performed, adopting 
the parametrized model of the prompt explosion from a collapsing 0-Ne-Mg core for the phys- 
ical conditions and compared with the results obtained with the HFB-7 and droplet-type mass 
formulas. Due to its weak shell effect at the neutron magic numbers in the neutron-rich region, 
the microscopic mass formulas (HFB-2 and HFB-7) give rise to a spread of the abundance dis- 
tribution in the vicinity of the r-process peaks [A — 130 and 195). While this effect resolves 
the large underproduction at A w 115 and 140 obtained with droplet-type mass formulas, large 
deviations compared to the solar pattern are found near the third r-process peak. It is shown that 
a solar-like r-process pattern can be obtained if the dynamical timescales of the outgoing mass 
trajectories are increased by a factor of about 2 — 3, or if the /3-decay rates are systematically 
increased by the same factor. 
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1. INTRODUCTION 

The origin of the rapid neutron-capture (r- 
process) nuclei is still a mystery. One of the 
underlying difficulties is that the astrophysical 
site (and consequently the astrophysical condi- 
tions) in which the r-process takes place has not 
been identified. Although some scenarios such 
as "neutrino-driven winds" from nascent neutron 
stars (Woosley et al. 1994; Takahashi et al. 1994; 
Qian & Woosley 1996; CardaU & Fuller 1997; Ot- 
suki et al. 2000; Wanajo et al. 2001; Thompson, 
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Burrows, & Meyer 2001), "prompt supernova ex- 
plosions" from small iron cores (Sumiyoshi et al. 
2001) or from 0-Ne-Mg cores (Wanajo et al. 2003), 
and "neutron star mergers" (Freiburghaus, Ross- 
wog, & Thielemann 1999) show some promise, 
each of them faces severe problems and cannot at 
the present time be called for to explain the pro- 
duction and galactic enrichment of the r-process 
nuclei observed in nature. 

Another underlying difficulty is due to the un- 
certainties in the theoretical predictions of nuclear 
data far from the /3-stability, for which essentially 
no experimental data exist. In particular, mass 
predictions for neutron-rich nuclei play a key role 
since they affect all the nuclear quantities of rele- 
vance in the r-process, namely the neutron cap- 
ture, photodisintegration and /3-decay rates, as 
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well as the fission probabilities . Although most of 
the recent mass formulas show fits to experimental 
masses of similar quality (characterized by an rms 
error of about 0.7 McV), the mass extrapolations 
far from the valley of /3-stability can differ from 
each other quite significantly (for a recent review, 
see Lunney, Pearson, & Thibault 2003). Recently, 
Hartree-Fock mass formulas with fully microscopic 
approaches have been constructed (Goriely, Ton- 
deur, & Pearson 2001; Samyn et al. 2002; Goriely 
et al. 2002; Samyn et al. 2003; Goriely et al. 2003). 
The latest Hartree-Fock-Bogoliubov formula, la- 
beled HFB-2 up to HFB-7(Goricly et al. 2002; 
Samyn et al. 2003; Goriely et al. 2003), are among 
the most accurate mass formulas, predicting the 
2135 measured masses with a root-mcan-square er- 
ror around 0.670 MeV for nuclei with iV, Z > 8, 
i.e with the same accuracy or even better than the 
one obtained with droplet-like mass formulas (Hilf 
et al. 1976; MoUer et al. 1995). A brief comparison 
of these mass models is given in § 2. 

The purpose of this study is to examine the 
effects of the newly-derived microscopic mass for- 
mula on the r-process nucleosynthesis and anal- 
yse to what extent a solar-like distribution can 
be obtained. Most of previous works devoted to 
the study of nuclear mass formulas and their im- 
pact on the r-process nucleosynthesis were based 
on site-independent approaches, the so-called 
canonical model assuming a constant temperature 
and neutron number density during the neutron- 
capture phase (Goriely & Arnould 1996; Cowan 
et al. 1999; Goriely & Clerbaux 1999; Goriely & 
Arnould 2001; Schatz et al. 2002). This approach 
is understandable since the astrophysical r-process 
site has not been unambiguously identified. An- 
other approach consists in considering promising 
sites and in modifying some of the relevant char- 
acteristics to force a successful r-process. In this 
case, the free parameter space is usually reduced 
and an easier analysis can be performed. In ad- 
dition, some works based on such site-specific ap- 
proaches have demonstrated that the abundance 
distributions are sensitive to the astrophysical con- 
ditions adopted, in particular during the "freeze- 
out" phase (Surman et al. 1997; Freiburghaus et 
al. 1999; Surman & Engel 2001; Wanajo et al. 
2002). 

We adopt, here, for the physical conditions the 
semi-realistic astrophysical model of the "prompt 



supernova explosion" from the collapsing 0-Ne- 
Mg core by Wanajo et al. (2003) (§ 3). The r- 
process nucleosynthesis with the HFB-2 mass for- 
mula in each outgoing mass trajectory is then cal- 
culated with a nuclear reaction network code. The 
mass-averaged yields over the mass shells relevant 
for the r-process is compared with the r-process 
abundance patterns in the solar system and in an 
extremely metal-poor star (CS 22892-052), as well 
as with those obtained with other mass formulas, 
more specifically the HFB-7 mass prediction and 
the extensively used droplet formulas of Hilf et al. 
(1976) and MoUer et al. (1995) (§ 4). Uncertain- 
ties stemming from the astrophysical conditions 
and extra nuclear ingredients are discussed in § 5. 
A summary follows in § 6. 

2. MICROSCOPIC MASS MODELS 

Among the ground state properties, the atomic 
mass is obviously the most fundamental quantity 
and influences the r-process abundance predic- 
tions mainly through the (n, 7) — (7, n) competi- 
tion taking place in the neutron-rich region. The 
calculation of the reaction rates also requires the 
knowledge of other ground state properties, such 
as the deformation, density distribution, single- 
particle level scheme, pairing force and shell cor- 
rection energies. Recently, impressive progress has 
been made experimentally. This situation results 
largely from recent measurements with Penning- 
trap or Schottky spectrometers which have en- 
larged the region of known masses, in particular 
towards the neutron-deficient side of the valley of 
nuclear stability (for a recent review, see Lunney, 
Pearson, & Thibault 2003). The new Atomic Mass 
Evaluation (Audi & Wapstra 2001) contains 2214 
measured masses, i.e 250 more than the 1995 one 
(Audi & Wapstra 1995). A more accurate mass 
determination is also available now for about 132 
nuclides originally included in the 1995 compila- 
tion. Out of the 382 new experimental masses, 
337 are located in the proton-rich region of the nu- 
clear chart and only 45 in the neutron rich region. 
As far as nuclei directly involved in the r-process 
are concerned, almost no experimental mass data 
exist, and theory must fill the gap. 

Attempts to estimate nuclear masses go back 
to the liquid-drop Weizsiicker mass formula. Im- 
provements to this original model have been 
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brought little by little, leading to the development 
of macroscopic-microscopic mass formulas, such 
as the droplet model (DM) (e.g Hilf et al. 1976) 
and the 'finite-range droplet model' (FRDM) of 
Moller et al. (1995). In this framework, the macro- 
scopic contribution to the masses and the micro- 
scopic corrections of phenomenological nature are 
treated independently, both parts being connected 
solely through a parameter fit to experimental 
masses. Despite the great empirical success of 
these formulas (e.g. FRDM fits the 2135 ^ > 8 
experimental masses (Audi & Wapstra 2001) with 
an rms error of 0.676 MeV), it suffers from major 
shortcomings, such as the incoherency of the link 
between the macroscopic part and the microscopic 
correction, the instability of the mass prediction 
to different parameter sets, or the instability of 
the shell corrections. As a consequence, its relia- 
bility when extrapolating far from experimentally 
known masses is severely limited. Although these 
models have been extensively used for astrophysics 
applications (Woosley et al. 1994; Takahashi et al. 
1994; Sumiyoshi et al. 2001; Wanajo et al. 2001, 
2002, 2003), there is an obvious need to develop a 
mass model that is more closely connected to the 
basic nuclear interaction properties. 

A new major progress has been achieved re- 
cently within the Hartree-Fock method (Goricly, 
Tondeur, k Pearson 2001; Samyn et al. 2002, 2003; 
Goriely et al. 2002, 2003). It is now demonstrated 
that this microscopic approach, making use of a 
Skyrme force fitted to essentially all the mass data, 
is not only feasible, but can successfully compete 
with the most accurate droplet-like formulas avail- 
able nowadays in the reproduction of measured 
masses. This holds true not only when the pair- 
ing force is described in the BCS approximation, 
but also when the Bogoliubov method is adopted 
(HFB model), which has the advantage of ensur- 
ing the self-consistency of the treatment of the nu- 
clear single-particle and pairing properties. These 
large-scale HFB calculations are based on the con- 
ventional 10-parametcr Skyrnic force along with a 
4-parameter (5-function pairing force. The Skyrme 
and pairing parameters are derived from a fit to 
the full data set of 2135 measured masses of nuclei 
with Z,N > 8, leading to an rms error of the order 
of 0.674 MeV for the HFB-2 mass table. 

Despite the success of the HFB-2 mass formula, 
a series of studies of possible modifications to the 



basic force model and to the method of calcula- 
tion was initiated all within the HFB framework 
in order to test the reliability of the mass extrap- 
olations (Samyn ct al. 2003; Goricly et al. 2003). 
For this reason, a set of additional 5 new mass 
tables, referred to as HFB-3 to HFB-7 were de- 
signed to analyse the sensitivity of the mass fit 
and extrapolations to the prescription used for the 
center-of-mass correction, the density-dependence 
of the pairing interaction and the amplitude of 
the effective nucleon mass. The new mass fits 
are of the same quality as the HFB-2 mass pre- 
dictions. In addition, it is found that globally 
the extrapolations out to the neutron-drip line 
of all these different HFB mass formulas are es- 
sentially equivalent. Figure 1 (lower panel) com- 
pares the HFB-2 and HFB-7 masses for all nuclei 
with 8 < Z < 110 lying between the proton and 
neutron driplines. Although HFB-2 and HFB-7 
are obtained with significantly different Skyrme 
forces (in particular, HFB-2 is characterized by an 
density-independent pairing force and an effective 
isoscalar mass M* = 1.04, while HFB-7 has an 
density-dependent pairing force and M* = 0.8), 
deviations smaller than about 2 MeV are obtained 
for nuclei with Z < 82. 

In contrast, more discrepancies are seen be- 
tween HFB and droplet-type masses (Figure 1), es- 
pecially for superheavy nuclei. For lighter species, 
the mass differences remain below some 5-10 MeV, 
although significantly different shell and deforma- 
tion effects are predicted. Most particularly, the 
HFB mass formulas show a weaker neutron-shell 
closure close to the neutron drip line with respect 
to droplet-like models as FRDM (for a more de- 
tailed discussion, see Goriely et al. 2002). Future 
improved microscopic predictions, in particular in 
mean field, but also shell model approaches, as 
well as future experiments will hopefully shed light 
on the behavior of the = 82 and 126 neutron 
shell closures at large isospins, in a similar way as 
for the A'' = 20, 28, and 50 magic numbers (Lun- 
ney, Pearson, & Thibault 2003). These major dif- 
ferences in the nuclear structure properties may 
affect the nucleosynthesis predictions, as shown in 
§4. 
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3. SUPERNOVA MODEL AND THE r- 
PROCESS 

Recent comprehensive spectroscopic studies 
of extremely metal-poor stars (Hill ct al. 2002; 
Cowan et al. 2002; Snedcn ct al. 2003), as well as 
chemical evolution studies (Ishimaru & Wanajo 
1999; Ishimaru et al. 2004; Argast ct al. 2003), 
in the Galactic halo suggest that the r-process is 
a primary nucleosynthesis process and that the 
astrophysics site could be associated with core- 
collapse supernovae. In fact, most of the recent 
studies on the nucleosynthesis r-process have been 
based on the "neutrino-driven wind" scenario, in 
which the free nucleons accelerated by the in- 
tense neutrino flux near the neutrino sphere of a 
core-collapse supernova assemble to heavier nu- 
clei (Woosley et al. 1994; Qian & Wooslcy 1996; 
Cardall & Fuller 1997; Otsuki et al. 2000; Wanajo 
et al. 2001; Thompson, Burrows, & Meyer 2001). 

In the present study, we use the parametrized 
model of the "prompt supernova explosion" from 
an 8 — IOM0 progenitor star (with a 1.38Mq O- 
Ne-Mg core) by Wanajo et al. (2003). The rea- 
son is that this model leads to r-abundance dis- 
tributions that have been shown to be relatively 
similar to the solar distribution, at least if an 
artificial enhancement of the shock-heating en- 
ergy is assumed. In addition, this scenario does 
not suffer from the problematic overproduction of 
A 90 nuclei seen in the neutrino-driven wind 
model (Woosley et al. 1994; Wanajo et al. 2001). 
Finally, the physical conditions in hydrodynami- 
cal (i.e., prompt) explosions may not be signifi- 
cantly affected by highly uncertain neutrino heat- 
ing (or cooling) processes (Hillebrandt, Nomoto, 
& Wolff 1984) as they are in the neutrino wind 
model. The site modeling is consequently easier 
and the final abundance predictions less sensitive 
to the neutrino physics. All these reasons justify 
the present study and the comparison of the nu- 
cleosynthesis results with the solar r-abundance 
distribution. It is obvious, however, that a simi- 
lar study based on the neutrino-driven wind sce- 
nario is needed since the prompt explosion sce- 
nario also faces severe problems that need to be 
resolved (Wanajo et al. 2003). In particular, it 
should be recalled here that in consistent realistic 
prompt explosion models, only a weak explosion 
is obtained and no r-processing. An artificial en- 



hancement of the shock energy is needed to drive 
an energetic explosion and provide the necessary 
conditions for a successful r-process. 

Wanajo et al. (2003) performed core-collapse 
simulations with a one-dimension, Newtonian hy- 
drodynamic code. The highly neutronized ejecta 
(YeOj initial electron mole fraction, w 0.14 — 0.20) 
in the (artificial) energetic explosion (model Q6) 
is subject to a strong production of r-process nu- 
clei. In this study, we use the mass trajectories 
of the ejected material from Mej = O.OSMq (zone 
number 60, Yeo = 0.45) to O.SIMq (zone num- 
ber 105, Yeo — 0.14) in mass coordinate (the sur- 
face of the 0-Ne-Mg core is at mass coordinate 
zero), as shown in Figure 2. This mass range was 
shown by Wanajo et al. (2003) to give a solar- 
like r-abundance distribution assuming that part 
of the ejected material with low- Ye (layers 106 to 
132 of Wanajo et al. 2003) was artificially assumed 
to be reaccreted onto the compact object, and con- 
sequently not to contribute to the galactic enrich- 
ment. The temperature and density histories of 
some mass elements are presented in Figure 3. 

Adopting model Q6 in Wanajo et al. (2003) 
for the physical conditions, the r-process abun- 
dances are obtained by solving an extensive nu- 
clear reaction network code. The network consists 
of ~ 4400 species, all the way from single neu- 
trons and protons up to the fermium {Z = 100) 
isotopes (Figure 4). We include all relevant reac- 
tions, i.e. (n,7), (^,7), (a, 7), (p,n), (a,n), (a,p), 
and their inverse. All reaction rates are calculated 
within the statistical model of Hauser-Feshbach 
making use of experimental masses (Audi & Wap- 
stra 2001) whenever available or the HFB-2 mass 
predictions (Goriely et al. 2002) otherwise. The 
photodisintegration rates are deduced from the re- 
verse (n, 7) rates applying the reciprocity theorem 
with the nuclear masses considered. The ,,9-dccay 
and /3-delayed neutron emission rates arc taken 
from the gross theory (GT2) of Tachibana, Ya- 
mada, & Yoshida (1990), obtained with the ETFSI 
(extended Thomas- Fermi plus Strutinsky integral, 
Aboussir et al. 1995) Q/j predictions. Other nu- 
clear inputs are the same as in Wanajo et al. 
(2003). 

A word of caution is required concerning some 
inconsistencies in the nuclear input considered 
here, and most particularly concerning the use of 
ETFSI masses to estimate /3-decay rates within 
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the GT2 model. As mentioned in § 2, nuclear 
masses influence the r-abundance distribution 
mainly through the (n, 7) — (7, n) competition tak- 
ing place in the noutron-rich region. Thanks to 
the balance theorem, such a competition remains 
relatively independent of the model adopted to 
estimate the neutron capture and photodisinte- 
gration rates. In contrast, the sensitivity of the 
absolute value of the reaction and /3-decay rates 
on nuclear masses can depend significantly on the 
model considered. For example, within the sta- 
tistical GT2 model, the /?-decay rate scale as the 
fifth power of the Qp value, while in the continuum 
QRPA model (e.g Borzov & Goriely 2003), the 
rates depend on the possible spin-isospin excita- 
tion within the Qp window but are not directly af- 
fected by the nuclear mass as such. Ultimately, for 
each mass model, a consistent estimate of the reac- 
tion and /3-dccay rates based on the corresponding 
mass model predictions for masses, but also de- 
formation, single-particle properties, . . . should be 
performed. In order to simplify here the discus- 
sion on the impact of masses, but also to avoid 
the complicate model dependence in the estimate 
of the reaction and /3-decay rates, we will restrict 
ourselves to consider in all cases only one set of 
/3-decay rates, namely the GT2 rates with ETFSI 
masses. 

The calculation in each mass trajectory is 
started at Tg = 9 (where Tg = T/10^ K). The ini- 
tial composition is taken to be that of the nuclear 
statistical equilibrium with the matter density p 
at Tg = 9, and consists mostly of free nucleons 
and alpha particles. The calculation is terminated 
after 10 seconds from the start in each mass tra- 
jectory. Snapshots of the r-process calculation in 
the mass-trajectory 99 (I.IIM0 in mass coordi- 
nate and YeO = 0.16) are shown in Figure 4. The 
abundances (mole fractions) are color coded, while 
the isotopes included in the reaction network are 
shown by dots with the stable and meta-stable 
isotopes by larger dots. Contours of constant 
two-neutron-scparation energies divided by two 
(S'2n/2) from 1 to 8 MeV are superimposed. The 
abundance curve as a function of the mass number 
is also shown in each upper left panel. 

The upper panel of Figure 4 shows the time slice 
when the neutron-to-seed ratio, Yn/Yf^, decreases 
to ~ 1, where Yn is the abundance of free neutrons 
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that of nuclei heavier than helium. At this time, 

the r-process flow closely follows an iso-52n curve 
given by S'2n/2 = 5° « 1.6 MeV (red hnc), where 
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Nn is the number density of free neutrons in units 
of cm~^. The condition Yn/Y^ = 1 is referred to 
as the neutron exhaustion. 

For the adopted (n, 7) rates, a quasi-equilibrium 
between neutron capture and photodisintegration 
continues until Tp/Tn decreases to ~ 1 as can be 
seen in the middle panel of Figure 4, where 
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are the averaged /3-decay and neutron-capture life- 
times, and X/3{Z, A) and A„(Z, A) the /3-decay and 
neutron-capture rates of the nucleus {Z,A), re- 
spectively. At this time, the abundance distribu- 
tion is still near the w 2.6 MeV line. We re- 
fer hereafter the condition Tg/r„ = 1 as freezeout, 
and the epoch between neutron exhaustion and 
freezeout as the freezeout phase. This freezeout 
phase is of particular relevance to determine the 
final abundance pattern as discussed in § 4. Af- 
ter the freezeout, the global abundance pattern is 
not significantly affected anymore, although it can 
still be somewhat smoothed out by photodisinte- 
gration and /3-delayed neutron emission, as can be 
seen in the lower panel of Figure 4 {t^/tu ~ 0.2). 

4. IMPACT OF MASS PREDICTIONS 
ON THE r-PROCESS 

The final r-process yields are mass-averaged 
over all mass trajectories between 1.08 and 
I.3OM0 in mass coordinate. In Figure 5a, the 
mass-averaged abundances (line) are compared 
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with the solar r-process abundance pattern (dots, 
Goriely 1999) that is scaled to match the height 
of the third r-process peak. For comparison, iden- 
tical calculations were performed by replacing our 
standard HFB-2 masses by the HFB-7 (Goriely et 
al. 2003, Figure 5b), FROM (MoUer et al. 1995, 
Figure 5c), and DM (Hilf et al. 1976, Figure 5d) 
predictions. 

Roughly speaking, an overall agreement be- 
tween calculated and solar abundance curves can 
be seen in Figures 5a-d for nuclei with A = 100 — 
200. Discrepancies common to all cases are, how- 
ever, observed, in particular the deficiencies for 
nuclei in the vicinity of the first r-process peak 
{A = 80), the shifts with respect to the solar curve 
of the second and third r-process peaks {A = 130 
and 195) to lower mass numbers 3 units in 
mass number) and the second to third peak height 
which is significantly smaller than observed in the 
solar distribution. We do not discuss further on 
the deficiencies near A = 80, since theoretically 
the results are highly dependent on the adopted 
mass cut of the ejecta (see § 3 and Wanajo et al. 
2003, for more details) and observationally there 
has been an increasing number of evidences that 
the light {A < 130) r-process nuclei have a differ- 
ent origin from the heavy {A > 130) ones (e.g.. 
Hill et al. 2002; Sneden et al. 2003). For the shifts 
of the peak locations, as well as the second to 
third peak ratio, it is conceivable that this defi- 
ciency does not originate from the mass model, 
but rather from inadequate astrophysical or nu- 
clear description, since it is observed for all mass 
formulas. This point will be further discussed in 
§ 5. 

A few significant differences in the abundance 
patterns can, however, be observed near the sec- 
ond and third peaks when use is made of the 
Hartree-Fock models (HFB-2 and HFB-7) on one 
side and the droplet models on the other side. 
First, the underproduction of nuclei at A w 115 
and 140 is more pronounced with the FRDM (Fig- 
ure 5c) and DM (Figure 5d) masses than with the 
HFB-2 (Figure 5a) or HFB-7 (Figure 5b) masses. 
Second, the abundances near A = 130 in the HFB 
cases are spread out in contrast to what is ob- 
served in the solar r-abundances. Third, the abun- 
dance curves near the third peak with the HFB 
masses are widened and the valley at ^ = 183 as 
observed in the solar r-distribution is significantly 



shifted to lower mass numbers. 

These differences reflect the model properties of 
iso-S'2n/2 curves shown in Figure 6, along which 
the r-process proceeds (see also Figure 4). Ma- 
jor local differences between the HFB (Figures 6a, 
b) and the droplet masses (Figures 6c, d) are 
found near the neutron magic numbers A'^ = 82 
and 126. The Hartree-Fock masses show weaker 
shell-closures, i.e. smoother iso-S'2ri/2 curves, at 
A = 82 and 126. This reduced shell effect (which 
actually is not totally quenched at the neutron 
dripline in contrast to other predictions, e.g Pear- 
son, Nayak, & Goriely 1996) is responsible for 
spreading the second and third abundance peaks. 
In contrast, the large droplet shell effect is clearly 
seen at A' = 82 and 126 by the stec^p 5'2n/2 char- 
acter even close to the neutron dripline and gives 
rise to sharp abundance peaks. 

Figure 7 illustrates the formation process of the 

second and third peaks with the HFB-2 masses 
in more details. The left panels show the snap- 
shots of the abundance distributions for the tra- 
jectory 90 (1.15M0 in mass coordinate and Fgo = 
0.23), and the right panels for the trajectory 99 
(I.IIM0 in mass coordinate and Yeo = 0.16) 
at the neutron exhaustion (top panels), freezeout 
(middle panels), and later time (bottom panels). 
The trajectory 90 (left panels in Figure 7) gives 
a major contribution to the second-peak forma- 
tion. The abundances are widely distributed al- 
ready at neutron exhaustion (left-top panel), ow- 
ing to the weaker shell gaps at A = 82, in par- 
ticular for low S'2n/2. This results in a split of 
the second abundance peak as seen in Figure 5a. 
The shift of the abundance distribution to the 
heavier mass number during the freezeout phase 
is small, although the abundance curve is sig- 
nificantly smoothed by photodisintegration (and 
slightly by /^-delayed emission). The reason is 
that the freezeout takes place at a low SjJ-value 
(« 2.9 MeV, left-middle panel in Figure 7). 

The formation process of the third peak in the 
trajectory 99 can be seen in the right panels of 
Figure 7. The A = 126 shell effect in the neutron- 
rich region is not reduced as the A^ = 82 one, 
though it still remains weaker than in the droplet 
approach. The splitting of the peak abundances is 
not observed, but the abundance curve is signifi- 
cantly broaden by photodisintegrations during the 
freezeout phase. The freezeout takes place when 
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is still low (2.6 McV), and for this reason the 
peak position at A = 190 is not shifted during the 
freezeout phase. 

Some of the deficiencies observed in the r- 
abundance distribution obtained with HFB masses 
are seen to be cured when considering droplet 
masses (Figure 5). In particular, the DM masses 
characterized by isospin-independent N = 82 and 
N = 126 shell effects give rise to an abundance 
curve which is seen to be globally in better agree- 
ment with the solar pattern than do the HFB 
curves. It would be tempting to judge the qual- 
ity of the nuclear mass predictions on the ba- 
sis of such a comparison. However, it must be 
recalled here that such r-process calculations in 
prompt explosion scenarios (like in any other site 
studied so far) remain strongly affected by other 
input uncertainties in the astrophysical model- 
ing as well as extra nuclear ingredients. More 
specifically, the total mass finally ejected (i.e not 
reaccreted onto the proto-neutron star) remains 
unknown and strongly dependent on the astro- 
physical modeling of the explosion. The same 
holds for the dynamical timescale of the outgoing 
material. Finally, major difficulties still exist in 
predicting reliably the /3-decay half-lives of exotic 
neutron-rich nuclei. These uncertainties are stud- 
ied in § 5, where it is shown that a modification 
of either the physical conditions of the explosion 
or the /3-decay rates can influence drastically the 
agreement with the solar r-distribution. 

Finally, it could also be argued that the out- 
come of the prompt explosion nucleosynthesis in 
one specific star might not need to match the so- 
lar r-process pattern. Recent spectroscopic stud- 
ies show that the elemental abundance patterns 
between the second and third peaks of r-process- 
enhanced ultra-metal-poor stars in the Galactic 
halo are in excellent agreement with the solar r- 
process abundances. In particular, the agreement 
extends to the third peak elements Os, Ir, and Ft 
for HD 115444, HD 126238 (Sneden et al. 1998), 
BD -M7°3248 (Cowan et al. 2002), and CS 22892- 
052 (Sneden et al. 2003, Figure 8), and CS 31082- 
001 (HiU et al. 2002). These observations led these 
authors to the conclusion that there is only one 
single r-process site that contributes to the pro- 
duction of heavy {A > 130) r-process nuclei and 
that the corresponding distribution is similar to 
the solar r-process pattern (Sneden et al. 1996). 



Instead of comparing the predicted abundances 
with the solar content, we can therefore use such 
observations and estimate to what extent the dis- 
agreements seen in Figure 5 and discussed above 
are also reflected in the comparison with the ob- 
served elemental distribution. The comparison is 
shown in Figure 8 for the most-studied, extremely 
metal-poor star CS 22892-052 (Sneden et al. 2003). 
Even in this case, it is difficult to discriminate be- 
tween the predictions obtained with the four dif- 
ferent mass formulas considered here. This is quite 
normal since the observed distribution is reduced 
to some 30 elements (in comparison with some 130 
nuclei in the solar material). And above all, the 
observed pattern could well be explained by the 
nuclear properties invariance in the 56 < Z < 76 
region rather than from specific astrophysics con- 
ditions as discussed in Goriely & Arnould (1997). 
It is therefore more informative to compare the r- 
abundance prediction with the solar pattern, keep- 
ing in mind that there is no stringent observa- 
tional constraint that the final distribution should 
be strictly solar. Nevertheless, future abundance 
determinations of as many extremely metal-poor 
stars as possible will be of particular importance 
to shed light on nucleosynthesis yields from only 
one or a few supernova events. 

5. UNCERTAINTIES IN THE r-PROCESS 
CALCULATIONS 

In this section, we examine possible astrophys- 
ical and nuclear modifications leading to a better 
agreement of the HFB-2 abundance curve with the 
solar r-process pattern. Regarding the prompt su- 
pernova explosions considered here, we do not ex- 
pect the entropy per baryon of the shocked ma- 
terial to be very diflferent from the one deduced 
from the simulation (~ 10 in units of the Boltz- 
mann constant), since many previous simulations 
resulted in similar values (Hillebrandt, Nomoto, & 
Wolff 1984; Sumiyoshi et al. 2001). However, the 
dynamical timescale of the outgoing matter (with- 
out changing the entropy) can significantly differ 
from model to model. For example, the dynami- 
cal timescales of the mass trajectories responsible 
for r-process in Sumiyoshi et al. (2001) are typ- 
ically a few times shorter than those applied in 
this study (Wanajo et al. 2003), because of their 
different treatments of electron capture. It is also 
conceivable that other effects, such as convection, 
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asymmetric explosion, or reverse shock from the 
outer envelope, cause an acceleration or a deceler- 
ation of the outgoing material. 

In order to test the impact of a change in the 
dynamical timescales, we simply modify the den- 
sity and temperature profiles of each trajectory, so 
that p'{t) = p{t/ft) and T'{t) = T{t/ft), i.e the 
dynamical timescale is multiplied by a factor of 
ft- Figures 9 and 10 show the abundance curves 
for ft ~ 0.2,0.5 (fast), 1 (standard), 2, 3, and 5 
(slow) obtained with the mass-trajectories 90 and 
99, which are responsible for the formation of the 
second and third peaks, respectively. The abun- 
dance curves near the second (Figure 9) and third 
(Figure 10) peaks are seen to agree fairly well with 
to the solar r-proccss pattern when the timescale 
is multiplied by a factor of 2 — 3 (Figures 9d-e and 
Figures lOd-e), and to disagree for the fast trajec- 
tories (Figures 9a-b and Figures lOa-b). 

This temporal effect is due to the significant 
abundance evolution during the feezeout phase for 
the slowly outgoing material, as shown in Fig- 
ure 11 (same as Figure 7, but for fi = 3). At 
the neutron exhaustion, the abundances are dis- 
persed near the second (left-top panel) and third 
(right-top panel) peaks, respectively, owing to the 
weak shell gaps at iV = 82 and 126 for low S'2n/2 
(« 2 — 3 MeV), as in the standard case (Figure 7). 
The freezeout takes place (middle panels), how- 
ever, at much higher values (4.2 and 3.7 MeV 
for the trajectories 90 and 99, respectively) than 
those in the standard case (2.9 and 2.6 MeV, Fig- 
ure 7). The reason is that the temperatures at the 
freezeout are higher (Tg = 1.4 and 1.2 for trajec- 
tories 90 and 99, respectively) owing to the slow 
expansion with ft — 3 than those in the standard 
case (Tg = 1.0 and 0.93). As a result, the abun- 
dance peaks shift to A = 130 (trajectory 90) and 
194 (trajectory 99), i.e mostly at the same posi- 
tion as those in the solar r-distribution. Concomi- 
tantly, the width of the abundance curve near the 
second and third peaks is reduced significantly due 
to the larger shell gaps at S'J^ w 4 MeV than that 
at S'JJ < 3 MeV, and becomes similar to the one 
observed in the solar pattern (Figures 9e and lOe). 

Assuming all the nuclear inputs are correct and 
the r-process indeed originates from the collapse of 
0-Ne-Mg core, this would imply that the expan- 
sion of the outgoing material must be significantly 
slower than obtained by Wanajo et al. (2003). A 



possible explanation could be found in the decel- 
eration of the ejecta by the reverse shock from the 
outer envelope if the hydrogen and helium layers 
survive the mass loss till the onset of the core col- 
lapse. This effect is absent from the simulation of 
Wanajo et al. (2003) with a bare 0-Ne-Mg white 
dwarf. This reverse shock effect would also reduce 
the total overproduction of r-processed material 
per event in the prompt explosion scenario (see 
Wanajo et al. 2003). It is interesting to note that 
such a significant deceleration of the outgoing ma- 
terial takes place in the neutrino-powered explo- 
sion due to the stalled shock wall once launched by 
a core bounce (Woosley et al. 1994). More detailed 
numerical simulations of the collapsing 0-Ne-Mg 
core including the outer envelope are needed to 
quantify this effect. 

The final mass-averaged abundance curve cor- 
responding to mass trajectories slowed by a factor 
of three {ft = 3) is shown with the scaled solar 
r-process abundances in Figure 12a. We find a 
good agreement between the calculated and solar 
r-process patterns, in particular near the second 
and third peaks. The steep valley at A ~ 183 
observed in the solar r-process abundances is also 
well reproduced. However, an underproduction at 
j4 » 115 and 140 appear. 

As far as nuclear uncertainties are concerned, 
/3-decays rates, neutron-captures rates and fission 
probabilities also play a key role (neutrino reac- 
tions are not important in the prompt explosion, 
see Wanajo et al. 2003). A change in their un- 
certain prediction can affect the nucleosynthesis 
predictions. In the present scenario, fission does 
not play an important role however, since it affects 
only a small fraction of the total r-processed mat- 
ter (Wanajo et al. 2003). The prediction of neu- 
tron capture rates, especially for exotic neutron- 
rich nuclei, remains difficult, in particular in the 
prediction of the direct capture contribution and 
of the 7-ray strength and neutron-nucleus opti- 
cal model potential for exotic neutron-rich nu- 
clei as required in the statistical reaction model 
of Hauser-Fcshbach (for more details, see Goricly 
1997, 2003). Much work remains to be done to 
estimate the neutron capture rates reliably. For 
this reason, their impact on the r-process nucle- 
osynthesis is postponed to a future study. 

Recent microscopic calculations of the /J-decay 
rate within the energy density functional plus 



8 



continuum QRPA formalism shows that the first 
forbidden transitions might have been underesti- 
mated in the past and can increase the /?-decay 
rate by a factor of about 2 along the neutron- 
rich N — 82 isotone and of about 3-10 along the 
neutron-rich A'' = 126 isotone compared with the 
calculation based exclusively on the allowed tran- 
sitions (Borzov 2003). Detailed microscopic esti- 
mates of the /3-decay rates (including forbidden 
transitions) for all nuclei of relevance in r-process 
calculations, as well as future experiments with 
radioactive ion beams are deeply awaited to im- 
prove this fundamental input in the description of 
the r-process nucleosynthesis. 

To estimate the influence of /3-decays, we show 
in Figure 12b the mass-averaged r-process yields 
obtained by multiplying all the /3-decay rates by 
a factor of three (i.e., reducing Tp by a factor of 
three). Interestingly, no significant differences are 
seen between Figures 12a and 12b. For /3-decay 
rates faster by a factor of three, the freezeout (cor- 
responding to T3 = T„) takes place at higher tem- 
peratures and thus at higher value, which has 
globally the same effect as slowing down the out- 
going material by the same factor. This test cal- 
culation shows the crucial role played by /3-decay 
rates when considering non-parametrized density 
and temperature profiles, i.e thermodynamic con- 
ditions fixed by (semi-)realistic models. 

6. SUMMARY 

We have examined the r-process nucleosynthe- 
sis with the state-of-the-art microscopic mass for- 
mula, HFB-2 (Goriely et al. 2002), to see its effects 
on the theoretical r-process abundance prediction. 
The yields of r-process species were calculated 
with a nuclear reaction network code (~ 4400 iso- 
topes), adopting favorable physical conditions in 
the prompt explosion model of a collapsing 0-Ne- 
Mg core by Wanajo et al. (2003). This model led 
to production of the solar-like r-abundance pat- 
tern in the previous study (Wanajo et al. 2003). 
However, an artificial enhancement of the shock- 
heating energy was needed to obtain requisite 
physical conditions for successful r-processing. 

Due to its weak shell effect at the neutron magic 
numbers in the neutron-rich region, the micro- 
scopic mass formulas (HFB-2 and HFB-7) give 
rise to a spread of the abundance distribution in 



the vicinity of the r-process peaks (A — 130 and 
195). While this effect resolves the large under- 
production at A « 115 and 140 obtained with 
droplet-type mass formulas, large deviations com- 
pared to the solar pattern are found near the third 
r-process peak. When use is made of droplet mass 
predictions, sharp r-process peaks are systemati- 
cally found, owing to their strong shell effect for 
neutron magic numbers even in the neutron-rich 
region. However, due to the numerous uncertain- 
ties still affecting the astrophysics models as well 
as the prediction of extra nuclear ingredients, it 
would be highly premature to judge the quality of 
the mass formula on the basis of such a compari- 
son. 

We found that abundance peaks similar to the 
one observed in the solar system could be recov- 
ered if the dynamical timescales of the mass tra- 
jectories are increased by a factor of ~ 2 — 3 (with- 
out any change in the entropy) or by decreasing 
systematically the /3-decay half-lives by the same 
factor. These changes might be conceivable when 
considering the current uncertainties in the astro- 
physics as well as in the nuclear /3-decay model. 
Much effort in the astrophysics and nuclear mod- 
eling remain to be devoted to improve the difficult 
description of the r-process nucleosynthesis. 
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for Scientific Research (13640245, 13740129) from 
the Ministry of Education, Culture, Sports, Sci- 
ence, and Technology of Japan. M.S. and S.G. 
are FNRS Research Fellow and Associate, respec- 
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Fig. 1. — Mass differences between HFB-2 and 
the DM of Hilf et al. (1976) (upper panel), HFB- 
2 and the FRDM of MoUcr ct al. (1995) (middle 
panel) and HFB-2 and HFB-7 (lower panel) as a 
function of the neutron number N for all nuclei 
with Z,N > 8 lying between the proton and the 
neutron driplines up to Z = 110 
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Fig. 2. — Distribution of the initial electron mole 
fraction in the ejected matter adopted in this 
study. The surface of the 0-Ne-Mg core is at mass 
coordinate zero. Selected mass elements are de- 
noted by zone numbers. 




Fig. 3. — Temperature (upper panel, Tg = 10^ K) 

and density (lower panel) histories for selected 
mass elements, which are denoted by zone num- 
bers. Time is set to zero at Tg = 9. 
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Fig. 4. — Abundance distributions at the neutron 
exhaustion (Yn/Yh ^ 1, top), freezeout {Tf^/Tn ~ 
1, middle), and later time (bottom) for the tra- 
jectory 99 (see text). The abundances are color- 
coded in the nuclide chart. The abundance curve 
as a function of mass number is shown in the upper 
left of each panel. The nuclei included in the reac- 
tion network are denoted by dots, with the stable 
and long-lived isotopes represented by large dots. 
Iso-S'2„/2 (= 1, 2, • • •, 7 MeV) and 5° curves are 14 
also superimposed. 
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Fig. 5. Final mass-averaged r-process abun- 
dances (line) as a function of mass number ob- 
tained with various mass formulas; (a) HFB-2, (b) 
HFB-7, (c) FRDM, and (d) DM. These are com- 
pared with the solar r-process abundances (points) 
from Goricly (1999), which are scaled to match the 
height of the third r-process peak. 
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Fig. 6. Contours of the S'2n/2 values (= 1, 2, • • •, 
8 MeV) for various mass formulas; (a) HFB-2, (b) 
HFB-7, (c) FRDM, and (d) DM. The stable and 
long-lived isotopes are also represented by dots. 
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Fig. 8. — Comparison of the mass- averaged yields 
(thick hnc) obtained with the HFB-2 masses, 
scaled at Eu [Z = 63), with the abundance pat- 
tern of CS 22892-052 (filled circles, with observa- 
tional error bars), as a function of atomic number. 
For Ga, Ge, Cd, Sn, and U, the observed upper 
limits are shown by the open circles. The scaled 
solar r-process pattern is shown by the thin line. 
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Fig. 9. — Final abundances as a function of mass 

nuinbcT from r-proccss calculations for the trajec- 
tory 90 with ft = (a) 0.2, (b) 0.5, (c) 1, (d) 2, (e) 3, 
and (f) 5 (see text). These are compared with the 
solar r-proccss abundances (points) from Goriely 
(1999), which are scaled to match the height of the 
second r-process peak. 
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Fig. 10. — Same as Figure 9, but for the trajectory 

99. The solar r-proccss abundances (points) from 
Goriely (1999) are scaled to match the height of 
the third r-process peak. 
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Fig. 12. — Same as Figure 5, but for (a) sfow tra- 
jectories {ft = 3) and (b) fast /3-decay rates (a 
factor of three, see text). 
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